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Abstract. The effects of radiative energy exchange on the convective 
instability of a weak field magnetic structure, which lead to a prediction 
and a physical explanation of the magnetic flux dependent field strength, 
are examined in detail using a real model stratification for the photo- 
spheric and convection zone structure of the Sun. Adopting the gener- 
alised Eddington approximation for the radiative transfer, which is valid 
both in the optically thick and the thin limits, we model the lateral ra- 
diative energy exchange by the tube with the external medium with a 
self-consistent inclusion of vertical radiative losses. 



1. Introduction 

One of the important properties of the small-scale solar magnetic structures, as 
established by the observations (Stenflo and Harvey 1985, Zayer et al 1989,1990; 
Schussler 1991), is that while the strong field network elements have field 
strengths very weakly dependent on the flux per element, 3>, with typical values 
of $ about 1 x 10 18 Mx and higher, the inner network weak field structures 
have a typical value of about 500 G(Keller et al 1994, Solanki et al 1996, Lin 
1995) for their strength and $ about 1 x 10 17 Mx or lower with a strong de- 
pendence between the two. If the convective collapse of a weak field tube is the 
global process responsible for the formation of the strong field tubes(Parker 1978; 
Webb and Roberts 1978; Spruit and Zweibel 1979; Spruit 1979; Hasan 1983, 
1984;Venkatakrishnan 1985; Steiner 1996), that comprise the network, with 
strengths weakly dependent on $ then it should be explained why tubes with 
smaller fluxes viz., the inner network elements do not collapse to kG strength; 
efficient radiative exchange with the surroundings by a small flux tube(Hasan 
1986; Venkatakrishnan 1986) offers a natural explanation; here, having in mind 
a quantitative comparison with the above mentioned observationally established 
properties of the Solar flux tubes, we employ a semi-empirical model of the pho- 
tospheric and the convection zone structure of the Sun and study in detail the 
effects of radiation on the convective instability and the wave motions. 



1 Joint Astronomy Programme, Indian Institute of Science, Bangalore 



1 



2. Equations 



We add to the familiar thin flux tube equations(Roberts and Webb, 1978) the 
following non-adiabatic energy equation(see e.g., Cox 1980), 



dp dp T\p 

dt dz p 



dp dp 
dt dz 



XT V.F (1) 



pc v T 



where p, p, and F denote the fluid pressure, density, and radiative energy flux 
respectively; (see Cox(1980) for the definitions of other thermodynamic quanti- 
ties); all the variables are evaluated on the tube axis(r=0); The radiative flux 
F is calculated in the generalised Eddington approximation following Unno and 
Spiegel(1967): the mean intensity J that is needed to evaluate the flux, 

Fr = -tt^VJ (2) 

6Kp 

is found by reducing the exact relation 

V.F R = Ak P {S - J) (3) 
to a form appropriate for a thin tube, which reads, 

Id 2 J 4/J e -J\ T „ 

- 3 -^ + i{— )= J - S <«> 

where dr = npdz, T a = npa, S is the source function which we take it to be 
the Planck function, a is the tube radius, and J e is the mean intensity in the 
external medium. J e is found by solving the equation 

1 ^ Je - T Q M 

2.1. Equilibrium 

The equilibrium stratification of the external medium that we use here is the 
one determined to match the combined VAL — C(Vernazza et al. 1981) and 
Spruit(1977) models (see Hasan,Kneer and Kalkofen(1998) for details about the 
construction of this external quite Sun model) ; the Rosseland mean opacities are 
calculated by interpolation from the tables of Kurucz(1993) for the upper layers 
and from those of Rogers and Iglesias(1992) for the deeper regions. We assume 
temperature equilibrium, T = T e , which implies that (3, defined as /3 = 8np/B 2 , 
is constant with z(if the dependence of p is neglected). The pressure and density 
are thus determined from 

V = — p Pe (6) 
P = — p Pe (7) 

The extent of the flux tube covers the atmosphere from the temperature min- 
imum in the chromosphere to 5000km deep in the convection zone with the 
photospheric surface(r = 1) being assigned z = 0. We measure the positive z 
downwards. 
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2.2. Linear Stability: The Perturbation Equations 

With the assumption that the perturbations in the ambient medium are negli- 
gible small perturbations inside the tube about the equilibrium described above 
obey the following equations in the linear limit: 
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(8) 
(9) 
(10) 
(11) 



where £ denotes the vertical displacement and H is the pressure scale height, 
the perturbation of equation (^) can be done in a straightforward manner: 



V.F 



dF, 



R 



R 



dz 



47TK aP (S' - J') + V.F R I^ + ^ 



(12) 



The perturbation in the mean intensity J is determined by perturbing and 
linearizing the transfer equation (||) which in the first order moment form are 




din J 7] din J J dlnJ TC 



dz rj 



dz J 



dz n 



(13) 




A f dln7i\n' A t a (1 + V c )(4 + 3t=?) j' 
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dz J n 
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qH„ T dz TC 



(14) 



where TC = F/4 is the Eddington flux, rj = up, H p is the pressure scale- 
height at the bottom and the various other quantities are as defined below: 



V, 



J 

s 



1 



is a measure of departure from radiative equilibrium, 

J — Jp 



(15) 



(16) 



is the ratio of the excess of the mean intensity inside the tube to the Planck 
function; e and q are the ratios, 



Tth 



(17) 
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Tth 

where 

m = ^ (19) 

is the radiative relaxation time over the length of one pressure scale-height at 
the bottom, 

_ pc^ _ IQaT 3 
Tr ~ K ' 3K p 1 Uj 

is the radiative relaxation time across the tube in the optically thick limit and 

™ = H < 21 > 

is the radiative relaxation time that one obtains in the optically thin limit (with 
Newton's law of cooling) (Spiegel, 1957); and, 

T a = K,pa (22) 

is the depth dependent optical thickness of the tube. We reduce the perturbation 
equations to a final set of four equations for the four variables £, p /p, J / J and 
Ti' /H. The optically thick reduction of the set of equations corresponds to 
taking the limit r tending to infinity and replacing the mean intensity by the 
Planck function. The optically thin case is got in the limit Tth approaching zero. 

2.3. Boundary Conditions 

We use closed mechanical boundary conditions and as thermal conditions we 
impose that there is no incoming radiation from above at the top boundary and 
that the perturbations are adiabatic at the bottom boundary. 



3. Results and Discussion 

The action of radiation that we bring out in this study of convective instability 
and wave-motions of the tube is explained conveniently with the help of the 
graphs shown in figs.{\) and (2) where we plot the growth rates and frequencies 
of the fundamental mode, which is the most unstable, respectively, as a function 
of the surface tube radius ao for various values of the plasma (3. 

3.1. Convective Instability 

The onset of the convective instability corresponds to the cusps in the curves 
of fig. (I) where the overstable mode's frequency becomes zero and the growth 
rate shoots up sharply. Comparison with the results obtained with only the 
lateral exchange in the diffusion approximation (Venkatakrishnan 1986, Hasan 
1986) or in the Newton's law of cooling reveal that these earlier treatments 
overestimate the degree of instability: the growth rates obtained with the present 
more accurate treatment of radiation with the inclusion of vertical exchange of 
radiation are appreciably smaller; moreover, for a given value of /3, i.e. for a 
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tube of given strength the onset of convective instability requires the size of the 
tube to be greater in the present case than that is required when the diffusion 
approximation is used. 

The convective instability is completely suppressed for tubes with the plasma 
(3 smaller than 2.45 whatever be its size; this corresponds to a field strength of 
about 1160Gauss at r = 1 inside the tube; this has to be compared with the 
value of 1350G(/3 = 1.83) that Spruit and Zweibel obtained in the adiabatic case. 
We point out here that the field strength of 1160G that we obtain here does not 
necessarily imply that all collapsing tubes of weaker fields will attain this unique 
value and become stable; this value represents a necessary strength for stability 
against convective collapse and thus can be considered as a minimum strength 
for stability; a collapsing weaker tube of sufficient size can of course attain an 
equilibrium collapsed state of field strength higher than this value(cf. Spruit 
1979). 

3.2. Size-strength Distribution for the Solar Tubes 

From the positions of the cusps in the curves of fig.(l), i.e. from the positions 
that mark the onset of the convective instability we pick up the values of the 
plasma f3 and the radius ao; the resulting radius- field strength dependence is 
shown in fig. (3); the corresponding flux-strength distribution in fig. (4). Com- 
parison of this curve with those observationally produced(Solanki et al 1997, 
Lin 1996) shows a remarkable agreement leading to the conclusion that indeed 
the convective collapse is the cause of the formation of the flux elements on 
the Sun's surface. Our refined, realistic, and exact treatment of the convective 
collapse process on the Sun reinforces the conclusions drawn from earlier sim- 
plified treatments (Venkatakrishnan 1986, Hasan 1986) and verifies the original 
suggestion and the physical explanation by Parker(1978) for the concentrated 
small scale magnetic structures on the Sun. 

3.3. Overstability 

The characteristics of the overstable mode are explained with the help of the 
fig.(l) again; the growth rates in the present Eddington approximation are 
lower than those obtained in the earlier treatments where the vertical losses are 
not taken care of and use either the diffusion approximation or the Newton's 
law of cooling; the differences thus demonstrate that overstability is hindered by 
the vertical losses 

We point out that while there is no damping out of oscillations when only 
lateral exchange takes place and the oscillations' growth rate only asymptotically 
goes to zero in the limit of large radii, i.e. in the adiabatic limit, the inclusion 
of vertical radiative losses make the oscillations to damp out for radii greater 
than a particular finite value which is determined by the /?; thus it is clarified 
that the horizontal exchange between vertically oscillating fluid elements acts 
to amplify the oscillations while the vertical losses always try to smooth out 
the fluctuations thereby introducing damping; this radiative damping is quite 
severe for the overstable mode of an intense flux tube on the Sun that it gets 
completely damped out for tubes of radii larger than a certain critical value. 

Finally we note that a tube which can undergo convective collapse for radii 
greater than a critical value for a given field strength remains overstable for all 
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Field strength at tau=1 inside the tube 



Figure 3. Size-strength relation for the stable tubes; see the text for 
details 
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Field strength at tau=1 inside the tube 

Figure 4. The magnetic flux-strength relation corresponding to the 
size-strength relation of the previous graph 
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smaller radii that it can take. 



4. Conclusions 

• We have demonstrated that radiative transport has a marked effect on 
the size-field strength relation for solar flux tubes. Our results can be 
be applied to solar flux tubes more reliably in view of the more refined 
treatment of radiative transfer. 

• We have generalized the necessary condition for the onset of convective 
instability in the presence of radiative energy exchange. We find that 
radiation has a stabilizing influence which is greater for tubes with small 
radius. 

• Overstability of the longitudinal slow mode is shown to be dependent on 
the tube radius: there is a critical tube radius above which a strong field 
convectively stable tube's oscillations get damped as a result of vertical 
radiative losses. 
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